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I. INTRODUCTION 



Many models of supersymmetry (SUSY) breaking predict that at least one of the SUSY 
partners of the top quark is sufficiently light such that pair production is possible at a future 
Linear Collider with cm. energies below 1 TeV. In such a scenario threshold studies are 
feasible, where squark pairs are produced with small relative velocities, and which allow for 
precise measurements of squark masses, lifetimes or couplings in close analogy to threshold 
measurements at the top-antitop pair threshold |2J, |3( . If the SUSY partner of the gluon 
is (as is general expected) not much lighter than the electroweak scale, then the low-energy 
QCD dynamics of squarks is, up to the fact that we are dealing with a colored spinless 
state, based on standard QCD with interactions carried by spin-1 gluons in the adjoint 
representation. 

It is a special feature of pair production of heavy colored particles close to the two-particle 
threshold that multi-gluon exchange leads to singular terms oc (a s /v) n and oc (a s lnv) n in 
the amplitude at n-loop perturbation theory, where v is the relative particle velocity. Thus 
for v ~ a s , which corresponds to a region in cm. energy of roughly \*/s — 2m\ < 0.1m, m be- 
ing the heavy particle mass, the singular terms have to be summed to all orders in a s . This 
is achieved most efficiently within an effective field theory. For top pair threshold produc- 
tion the effective theory vNRQCD ("velocity non-relativistic QCD") (HIHEf was developed 
to carry out this summation program aiming for NNLL order precision, i.e. accounting for 
0{a 2 s) a s v, v 2 ) qcd 0, I, 3 as well as 0(a em ) electroweak corrections For squark- 
antisquark pair threshold production previous studies involved leading-logarithmic (LL) 
precision 0, 0|. Higher order QCD or electroweak effects were not taken into account 
in a systematic manner, although it is known from the top threshold that corrections to the 
LL approximation can be substantial. 

In this work we present the vNRQCD Lagrangian for a particle-antiparticle pair of 
heavy non-relativistic colored scalars in the fundamental representation of QCD. We provide 
matching conditions and anomalous dimensions required for renormalization group improved 
computations of threshold pair production at next-to-leading logarithmic (NLL) order and 
scalar-antiscalar bound state energies at NNLL order in the non-relativistic expansion. For 
our presentation we follow closely the conventions and notations of Refs. 0, E|. Among the 
crucial elements of the construction are that vNRQCD is obtained by a single matching 
procedure at the hard scale m and that the soft and the ultrasoft renormalization scales, 
fis and /i;/, are correlated according to the non-relativistic energy-momentum relation of a 
heavy particle at all times. Thus we have = [Ag/m = mu 2 , where v is the dimensionless 
renormalization group scaling parameter of vNRQCD. (For corresponding computations in 
the pNRQCD approach |13( see Ref. Q.) For convenience, for the most part of the paper, 
we call the scalars in the fundamental representation of QCD simply "squarks" and the 
scalar versions of vNRQCD and QCD frequently just vNRQCD and QCD, respectively. 

The outline of the paper is as follows: In Sec. |H] we present the terms of the scalar 
vNRQCD Lagrangian relevant for this work together with their matching conditions. A 
special set of operators generated by ultrasoft renormalization is discussed in Sec. ITTL>. The 
running of all relevant potentials is discussed in Sec. IIII1 Finally, Sec. |IV| is devoted to 
the anomalous dimensions of the S- and P-wave squark-antisquark production currents. 
Conclusions are given in Sec. 
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FIG. 1: Effective theory diagrams for the interaction of a heavy squark with an ultrasoft gluon. 



II. EFFECTIVE THEORY LAGRANGIAN 

The scalar vNRQCD effective Lagrangian is written in terms of the fields ^ p for the 
squark, Xp for the antisquark, (A%, c q , (p q , <f> q ) for soft (massless) gluons, ghosts, quarks and 
squarks and (A^ 1 , c, (p, 0) for ultrasoft (massless) gluons, ghosts, quarks and squarks. For 
the light quarks and squarks we assume rif and n s light flavor components, respectively. 
The covariant derivative is = + ig u A^ = (D°, -D), with D° = d° + ig u A°, D = 
V — ig u A and only involves the ultrasoft gluon field and the ultrasoft gauge coupling. For 
the representation of the (anti)squark fields and the soft fields the vNRQCD label formalism 
is used to separate the fluctuations at the scales mv (soft) and mv 2 (ultrasoft) 0. Thus 
the spatial dependence of fields refers only to ultrasoft fluctuations and, for example, for 
the heavy squark field we have d^ippix) ~ mv 2 ip p (x), while the label p refers to the soft 
three-momentum component of the squark field. 



A. Basic Lagrangian 

The basic vNRQCD Lagrangian can be obtained by tree level matching to full QCD and 
can be separated into ultrasoft, soft and potential components, C = C u + C s + C p } The 
ultrasoft piece of the effective Lagrangian describes the interactions of ultrasoft gluons with 
squarks (Fig. and also contains the kinetic energy contributions for the squarks and the 
ultrasoft fields. It has the form 



iD 



o ( P -*D) 2 , p 4 



2m 



+ ^"3 + ■■■ 



i=l 



i=l 



(1) 



where G£ u is the ultrasoft field strength tensor and m is the heavy squark pole mass. The 
spatial dependence and the summation over the color indices of the fields are understood 
here and for the rest of this paper. To avoid double counting the interactions have to 
be multipole-expanded with the scaling p ~ mv and D ~ mv 2 . The form of the first 
term of Eq. (0) is determined by reparametrization invariance 3> IB an d does not receive 
any non-trivial renormalization factors. For the squark fields we use the non-relativistic 



normalization with ipr 



m 



3/2 



such that the squark kinetic energy has the usual form 



known from NRQCD for heavy quarks and non-relativistic quantum mechanics. To the order 



1 We suppress throughout this paper the renormalization constants that relate bare and renormalized quan- 
tities. Note that we use the convention that the antiquark field describes positive energy antisquarks. 
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FIG. 2: The Compton scattering diagrams a-g in the full theory generate the effective theory soft 
gluon coupling h. The diagrams e, f and g describe interactions with soft quarks and squarks and 
ghosts, respectively. Soft gluons, ghosts, quarks and squarks in the effective theory are collectively 
depicted by zigzagged lines. 

we are working C u agrees with the ultrasoft Lagrangian for the quark case. Likewise, the 
momentum and v scaling of the effective theory squark fields and effective theory Feynman 
diagrams is equivalent for the heavy quark case discussed in detail in Refs. 0,0]. Here we 
employ the same power counting formalism. The factors of fifj and //| that appear in the 
effective Lagrangian in d = 4 — 2e dimensions are uniquely determined by the requirement 
that the kinetic terms in the vNRQCD action are of order v 0, 0|. To be specific, for 
the squark and gluons fields this gives the scaling ip p ~ (mt>) 3//2_£ , ~ (mu 2 ) 1-6 , and 
A't ~ (m«) 1_f . For the covariant derivative this leads to the form = <9 M + i\L\ J g u A i> ' in 
dimensional regularization. 

The soft Lagrangian describes interactions of soft (massless) gluons, ghosts, quarks or 
squarks with the heavy squarks and has terms 



£ s = { - \°s" G % + c q q 2 c q + J2 i vm + E C 

q i=l i=l 

[Cq>, Cq}Y^ + J2Wp> T BZ ^ V>p) {^qn'T B ^q) 



(2) 



1=1 



1=1 

where Gf is the soft gluon field strength tensor and g s = g s {ns) (l^s = mv ) is the soft QCD 
coupling. We also allow for n s massless scalar flavor fields <p q . The theory graph describing 
the soft interactions is shown in Fig. |2h- One- loop time-ordered products of two soft vertices 
with a and a' in a four-squark matrix element contribute to scattering diagrams at order 
v a+a '~ l . The tensors uf%) , wffi, Y^ a \ and depend on the labels p',p, q,q' and 

can be obtained from matching to the full theory diagrams in Fig. |2K-g- Their explicit form 
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in Feynman gauge and only employing the soft gluon momentum label q is 

„<o) 1 „(0) _(V-2p-qy ( o) (p - p' - q)' (0 ) -2g°<fo' 
^oo - g0 , v 0t - - (p/ _ p)2 , ^ - (p/ _ p)2 , u tj - (p/ _ p)2 , 



(p' + p)-q (p' + p)-q , (P' 2 -P 2 ^ 



r/W _ _ _ _ ^ 

00 '"' 2m{q ) 2 m(p' - p) 2 ^ 2m(p' - p) 2 ' 

„(D (p+pT , ?°(p+p'r 



0i "~ 2mg° 2m(p'-p) 2 ' 

„d) = _ (p+pt , g°(p+p ; r 

10 2mg° 2m(p' - p) 2 ' 

E#> = [2 < PV + ^ m (2p / -2p-q) J ' + y m (p-p / -q) < ] 
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r/ (2) = CgQy - p) 2 (p ■ q) 2 + (p' ■ q) 2 (1 ~ c D )q 

00 8m V 2m 2 (g ) 3 Am 2 q° 

| (2-c D )(p-p')-q 

4m 2 g° 

„(2) [p-q(2p + q) t + p^q(2p / -q) t ] (c D - l)(p - p' + q)' 



4m 2 (g ) 2 4m 2 
c D (2p'-2p-q)* (P"-P 2 ) ( ^ 
8m 2 4m 2 (p'-p) 2lP P; 



(2p'-2p- q y(p' 2 -p 2 ) 2 
4m 2 (p' -p) 4 

rr (2) _ [p ■ q (p + p' + q)' + p' ■ q (p + p' - q)' ] (cd - i)q* 

<AT> r,/ ONO "T 



io - 4m 2 (g ) 2 4m 2 

cj(p - p' - q) j _ (p' 2 ~p 2 ) , ,y 
8m 2 2m 2 (p'-p) 2lP P) 

(p-p'-gHp' 2 -? 2 ) 2 
4m 2 (p' - p) 4 

„(2) = (p + p')'(p + p') J ' q'fr-p' + q)' _ g°^(p /2 - P 2 ) 2 
H 4m 2 g o + 4m 2 g o 2m 2 (p' -p) 4 ' 

^ = o, 

w (i) _ J_ (p - pO • q ^d) _ (p - p' + q)' w (d _ .zSl _ ^ /3) 

W °° ~ 2m + 2m(g0) 2 ' ^ " 2mg° ' W * ~ 2mq° ' ^ " 2m {6) 
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y(0) = ~g° y(l) = q • (P + PQ 

(p'-p) 2 ' 2m(p'-p) 2 ' 



y(2) 



c D q° 
8m 2 



z (o) = 1 z (0) = Z (i) _ n _ ~(P + P / Vt 



° ~ _ (P'-P) 2 ' ~2m(p / -p) 2 ' 

7 (2) _ C D (2 ) _ 

X£ ] = Z<£\ a = 0,l,2. (4) 

Note that for completeness we also display terms proportional to p' 2 — p 2 . We emphasize 
that we do not use these terms in the definition of the soft Feynman rules. In time-ordered 
products of two soft vertices these terms vanish for on-shell scattering amplitudes. As can 
be checked from Eq. (J2J), the tensors involving soft gluons satisfy the relations 

u»M) = -UvM) > w,M) = w v ,{ q ') > q' = -q+p'-p, 

for on-shell external momenta and up to terms proportional to (p /2 — p 2 ). The expressions 
for o" = agree exactly with the corresponding results for the heavy quark case 0,0]) an d 
the expressions for a = 1 also agree up to the spin- dependent contributions that are absent 
for scalars. This is related to the fact that the spin independent part of the 0(1/ m) HQET 
(heavy quark effective theory) Lagrangian agrees with the effective Lagrangian for a single 
heavy squark, for simplicity called HSET (heavy scalar effective theory) in the following. 
This is because the spin independent parts only contain kinetic energy contributions which 
are fixed by reparametrization invariance 0, [l6[ and do not pick up non-trivial renormal- 
ization factors. Since ultrasoft corrections are only accounted for in the renormalization of 
time-ordered products of at least two soft vertices describing squark-antisquark scattering, 
the coefficients in C s are only renormalized due to soft interactions 0,^^. In analogy to the 
quark case the structure of the resulting soft divergences is identical to the one in HSET. 
Thus one can determine the leading-logarithmic (LL) soft running of the coefficients in C s 
by scaling the HSET Lagrangian at order 1/m 2 , 

£hset = r{iD° + ^ + c D g £iS + Q {—^ U , (5) 

to /i = mv, and then by matching the soft vNRQCD vertices in C s to the corresponding 
HSET diagrams. This leads to the dependence on g s and Co shown in Eqs. (J2J and 0- Note 
that in Eq. JSJ) the covariant derivatives refer to soft gluon interactions. For the solution for 
the coefficient cd{v) for heavy scalars we find 



/20 32CV 

Cd{v) = U + ucl 

with 

a s (m) 



z -13C A /(6Po) 



(6) 



z = 



fio= 1 jc A - i -^T F (7) 
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and the group theoretical factors Ca = 3, Cf = 4/3, Tp = 1/2. The term fi is the one- 
loop QCD beta-function that also accounts for the running due to the light squark flavors. 
Note that the coefficient cd in HSET vanishes at the hard scale, cd(1) = as shown in the 
appendix 1X1 where the tree level HSET Lagrangian is derived. The evolution of cd can be 



deduced from the results given in Ref. 3, 19j and differs from the corresponding coefficient 
in HQET because for the heavy scalar field the mixing to spin-dependent operators does not 
occur. The relation of the tensors in Eqs. to the HSET Lagrangian also explains (except 
for the dependence on the spin-dependent HQET coefficients cp and cs which are absent in 
the scalar case) the close resemblance of the a = 2 tensors to the quark results obtained in 
Ref. 0, 0] . We note that in Eqs. all ie prescriptions in the l/q° terms originating from 
integrating out the intermediate static squark are dropped. This prevents pinch-singularities 
from the l/q° term in one-loop diagrams related to time-ordered products of soft vertices. 

The potential Lagrangian describes potential-like squark-antisquark scattering interac- 
tions and has terms 



p>p' 

+ £ t4% ABC (p, p) Mi*? ) [r P .T A Tp pX t p ,T \- p 
p,p' 



where 



V{p, p') = (T A ® T 



vP | v^v | y^(p 2 + p' 2 ) | 



+ 



,(T) 



(8) 



-,ABC 



(P)P') 



m|k| 

;W^2 



2m 2 k 2 



1 1 h . . . 

k 2 mk rn? 



k 1 



ABC 



(9) 



and k = p' — p. The terms with the coefficients V c 1 are the Coulomb potential that 

( 1 T~) 

contributes at order a s v _1 to scattering diagrams. The terms with the coefficients V2 ' 

(T) 

and Vr contribute at order a s v and are the scalar analog of the Breit-Fermi potentials 
known from QED. The terms with the coefficients ' are only generated by one- loop 
diagrams and contribute at order a 2 s v°. For scattering of squark-antisquark pairs in 
the color singlet state only the linear combination of coefficients v\ = —CpV^ + V$ 
is relevant. Graphically the first and second terms in Eq. (JSJ) are depicted in Fig Ell and 
e, respectively. Matching to the full theory Born diagrams in Fig. |3k,,b leads to the LL 
matching conditions at v — 1, 



Vf >(1) = 4vra s (m) , Vf \l) = 4vra s (m) , Vf = -na s {m) 

v«(i) = o, v?\i) = o. 



(10) 



These matching conditions are consistent with the results given earlier in Ref. [20j. Note 

(T) 

that the matching condition for V2 differs from the corresponding result in the quark 
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FIG. 3: (a-c) QCD tree level diagrams relevant for the LL matching conditions of the potentials; 
(d,e) effective theory Feynman diagrams for the potential interactions. 



case. The annihilation diagram in Fig. does not contribute to scattering for a squark- 
antisquark pair in a color singlet state. Moreover, the annihilation diagram contributes at 
order a s v 3 because a squark-antisquark pair that annihilates into a gluon is in a P-wave 
state. Although the results will not be needed at the order we are interested in, we give 
them here for completeness: 



k(p,p') = (r A ®r 



VP p'.p 



+ . . 



rrr 



+ (1 



V« p'.p 



+ . . 



rrr 



where 



VP 



c 



-vra s (m), VP 



9i 



7ca s (m) 



(12) 



(1 T) 

The matching conditions for the coefficients ' arise from the terms in the one-loop full 
theory scattering diagrams contributing at order a 2 v°. Since in the effective theory all one- 
loop contributions at order a 2 v° from time-ordered products of soft vertices with o = 
and a' = 1 vanish the matching conditions are just equal to the corresponding full theory 
contributions. The latter can be obtained using the threshold expansion and read 



V, 



a 2 (m) 



7Ca C a 



Vl 1 \l) = ai(m) 



where 



Cd — 8Cf — 3Ca 



Ci 



(13) 



(14) 



The results agree with the corresponding matching conditions for the quark case 0, 
Similar to the agreement for the soft a = 0, 1 vertices this feature can again be traced back 
to the equivalence of the HQET and HSET effective actions mentioned above. 



B. Operators Generated by Ultrasoft Renormalization 

In the potential and soft components of the vNRQCD Lagrangian there are a number of 
operators with Wilson coefficients that have vanishing matching conditions at v = 1, but 
which become non-zero for v < 1 due to ultrasoft renormalization 0. Since the structure of 
these operators only depends on the form of the ultrasoft Lagrangian C u and on operators 
in the soft and potential Lag r&ngiciiis >^s,p t licit can contribute to scattering terms at order 
v~ x , they are in complete analogy to the corresponding operators in the quark case already 
discussed in detail in Ref. apart from the operators involving soft squarks which were not 



considered there. For completeness we briefly discuss the full set of these operators relevant 
to the order we are interested in. 

In the soft sector we have to include 6-field operators describing squark-antisquark scat- 
tering with two additional soft massless gluons, ghosts, quarks or squarks. The additional 
terms in the soft Lagrangian have the form 



r (2) n (2),(l) , r (2) n (2),(T) r (2) n (2),(T) 



J 2b ^2b 



J 2c ^2c 



(15) 



where 



r o(2),(i) 

^2a 



(0),(1) 
2i 



o 



(2),(T) 
2b 



m 2 2^ 



(0),(T) 



o, 



(2),(T) 
2c 



(P 2 + P /2 ) 



(T) 



(16) 



and the sums are related to the operators 



O 



(o) 

2 r - 
,(0) 



<? s 4 /4 e (V£ C ^p) (x-p' r Sx x-p) rw ? 



9/ ' 



= gUs (€> rZ V> P ) (x-p' x- P ) (^r^^) 



^ = <7 4 /4 e (V£ r^^p) (x V n u i x-p) (0 9 r 

,4 ..4e 



(0) 



,(0) 



'q) i 



<0) 



(0) 



,(0) 



0. 



(0) 
2c 



(V£ rgj Vp) (x* P > rg x-p) (c, i?> 



(17) 



referring to soft quarks, squarks, gluons and ghosts, respectively. Note that in Eq. (fTKj) 

(2) 

we have suppressed all sums over soft labels and flavor indices. The operators 2i are 
represented graphically by the diagram in Fig. H^. The matrices Tj are functions of the 
heavy squark and soft labels and have the form 



r(o),(i) 



,(0),(T) 



-(OUT) 
<j>,ij> 

<0),(T) 



<0),(T) 

A,i> 



F (0),(T) 
1 c,V 



F (0),(T) 
1 c,V> 



F (o),(T) = i 
* 2 



(2g° 7 ° + k ■ 7) (2g° 7 ° - k ■ 7) 



k 2 + 2k ■ q 



rpA 

F (0),(T) 

r 



p(0),(T) 



^(O).(l) 



k 2 - 2k • q 

1 r(o).(i) . 

' v,x 



r(0) 



+ 



{4 



0N2 



k 2 + 2k • q k 2 - 2k • q 



F (0),(T) 

(o),(i) = Ca 

A,fj,is r\ 



.(0),(1) 



^(o),(i) 



(0) 



£>#?(-?, p,p') (^ (0) )%(g,- P ,-pO 



1 



k 2 + 2k ■ q 



+ 



tf$(?,P,p') (f/ (0) )%(-g,- P ,-p') 



k 2 - 2k ■ q 



T- 



r (o),(i) 



^0),(T) 
= ^4 



1 



^(o),(i) 
A,i> 



+ 



F(0),(T) 



k 2 + 2kq k 2 

_ T^A r (0),(l) _ 

- 1 > 1 c,V ~ 



2k - q 



F(0),(1) 



(18) 
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FIG. 4: (a-f) Ultrasoft diagrams with p. A vertices for the computation of the LL anomalous 

(2) (1 T) 

dimension of the operators 2i ' (wave function renormalization graphs are understood); (g) 

(2) (i J 1 ) 

graphical representation of operators 2i ' , where i = a,b, c. 



The color indices for the two soft particles have been contracted in the expressions in Eq. (|18|) 
because at the order we are interested in we only need the operators in terms of one-loop 
diagrams where the soft lines are closed up (see Fig. [BJi) . When the soft lines are closed up 
the operators contribute to scattering matrix elements at order a s v. The Wilson coefficients 
run due to UV-divergences from ultrasoft gluon loops (and the corresponding soft pull-up 
terms) dressing the time-ordered products of two soft vertices with a — a' — 0, see Fig. HI 
The results have been worked out in Ref. 16] and read 



n {2), \ 4Ci ( 2 ). 3Ca -C d -4:C F 

6po ofJo 

<£\v) = InH , (19) 



where 



w = as{mv2 } . (20) 

In the potential sector we have to account for additional 4-squark operators which in- 
volve sums over intermediate 3-momentum labels. The additional terms in the potential 
Lagrangian are 

AC - V {1) {1) + V {T) (T) + V (1) (1) + V (T) (T) + V {T) {T) (21) 
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FIG. 6: Diagrams for the computation of the LL anomalous dimension of the order a s v poten- 
tials: a-f) ultrasoft gluon graphs with p. A vertices (wave function renormalization graphs are 



understood), g) soft graphs with a + a' = 2, h) soft graphs involving the operators O 



(2) 
2i 



before the remaining sums over the squark-antisquark soft 3-momenta labels are carried 
out. The operators O ki ' arise from the order a^v diagrams Figs. EH-c and are graphically 
depicted by Fig. EJ. The operators 0^' T ^ arise from the order a\v _1 diagrams Figs. EtLe 
and are graphically depicted by Fig. |I)g. The Wilson coefficients were computed in Ref. p 
and read 



v£V) 



v<TV) 



8CaC± 



ln(w) 



3/3 

CaC\{Ca + Cd) 



8C A 
Wo 



3/3 
Ci + 



V, 



(T). 
c2 1 



V) 



2Ca{Ca + Cd 
Wo 
16C A C 



ln(u>) 



3Po 



1 w ^ 
- in{w) 



(C A + C d ) 
32 



ln(w) . 



(24) 



Explicit expressions for 4-squark matrix elements of the operators and Od with the sum 
over labels carried out in n = 3 — 2e dimensions can be found in the appendix of Ref. . 

Note that it is possible to choose a scheme where the 14 potentials in Eq. © are also 
written as sum operators in analogy to Eq. (J2~T| . In this scheme the UV- finite contributions 
in NNLL matrix elements and the UV- divergences in NNNLL matrix elements differ. More 
details on this scheme are given in Ref. 0. 



III. ANOMALOUS DIMENSIONS FOR POTENTIALS 

The effective theory one-loop graphs required to determine the LL anomalous dimension 
of the order v potentials in Eq. (jUJ) are shown in Fig. |HJ Note that for ultrasoft gluon 
diagrams only the p. A couplings have to be accounted for because the contributions from 
the time-like A gluons vanish [17]. This feature can be formally understood from the fact 
that the leading order coupling of the heavy particle with time-like gluons can be removed 
by field redefinition |2lj . We find the anomalous dimensions 
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8 32 
-2(A) + -C A ) af(mi/) + yC A « s Ha s (mi/ 2 ) - 4 / 3 cf ) (z/) a£(mi/) 



A) 

2 
4 



1 + c D (i/) 



+ 9a 

6 



28 - lWy 



JO, 
6 



+ ^(4C F + C d - 3C A ) a s (mis)a s (mv 2 ) - 2ft Cg>) a 2 (mu) 



— Ci a 2 {mv) - y C x a s {mv)a s {mv 2 ) -2f3 C^J{u) a 2 {mu) 



a s {mv) 

r(2) ( 
2b 1 

t(2). 



(25) 



The results for and agree with the results for heavy quarks obtained in (and are 
consistent with corresponding results in Ref. [H]), while the result for differs due to 
the absence of the spin dependent HQET Wilson coefficient cp. Concerning the structure of 
logarithms at order a 2 described by Eq. fl23j) we disagree with the results given in Ref. jiO ]. 
Note that part of the discrepancy is related to approximations made in Ref. that are 
unjustified in the NRQCD v power counting scheme. (See also Ref. |H for more details.) 
Accounting for the matching conditions shown in Eqs. (jl(Jj) the solutions for the coefficients 
read 



Vf>( 



V. 



en, 



+ 



+ 



47r a s (m) z - 32? ^ A a s {m) z \n{w) , 
3p 

— ttz a s (m) 

7T [C A (152C F - 35C A ) - 4/3b(5CU + SC F ) 

13PoC A 
8tt(3/3o- 11C a )(5C a + 8C f ) 
13C A (6P - 13C A ) 

4tt(4CV - 2C A 



(26) 



■aJm) (z - 1) 



a s {m) [ z i-uc A /m) _ X ] 



4ttCi 



a s (m) (1 - z) + 



ot s (m) z \n(w] 
167rCi 



3{3 



a s (m) z ln(u?) . 



(27) 



(1 T) 

The running of the coefficients is determined solely by soft diagrams and associated 

with the known running of the strong coupling |6j . The results for the anomalous dimensions 
relevant for describing the squark-antisquark dynamics at NNLL order are 



-2 



^(^) + A^ + /5,^ (mi/) 



47T 



(4tt) 



~iyC A -C d ) [a s {mv)f 

Z7T 



where /?o,i,2 are the coefficients of the QCD beta function (in the MS scheme for ft). 
The anomalous dimension for the Coulomb potential is needed at three loops because the 
Coulomb potential already contributes at order a s t> -1 . The results agree with the results 
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(a) V r _ 2 (b) V kU2 (c) V k 




FIG. 7: Three-loop current correlator diagrams and the counterterm diagram for the computation 
of the NLL anomalous dimensions of cs,p- 

obtained in the quark case For the running of the Coulomb potential the agreement 
is obvious since the respective leading order interactions of heavy squarks and quarks with 
low energy gluons are universal. The agreement for the order v° potential coefficients V J [ 1 ' T ' 1 
is again related to the equivalence of the HQET and HSET effective actions at order 1/m 
discussed before. With the matching conditions shown in Eqs. (IT0|) and (fT3|) the solutions 
read 

VP{U)= {7CA - Cd) c&mv), V?V) = y c? a {mu) , Vf > {v) = Anaf (mv) , (28) 
where af^mu) is the QCD coupling with 3-loop running. 



IV. ANOMALOUS DIMENSIONS FOR PRODUCTION CURRENTS 

We consider the currents describing squark-antisquark production in S- and P-wave states 
relevant for 77 and e + e~ collisions, respectively: 

Js, P =r P x- p , j P , P = r P Px- P , (29) 

where cg(z/) and cp(v) are the corresponding Wilson coefficients. Concerning the summation 
of large QCD logarithms the LL matching conditions are conventionally normalized to unity. 
For e + e~ annihilation the NLL order matching condition for the P-wave vector current can 
be obtained from matching to the one-loop QCD amplitudes. From the QCD one-loop 
results obtained in |23j, |24| we find 

cp(1) = 1 - CFas{m) ■ (30) 

7T 

The corresponding matching condition for the S-wave current relevant for 77 collisions is 
presently unknown. In analogy to the quark case the currents in Eqs. (I29|) do not run at LL 
order since there are in general no UV-divergences at one-loop in the effective theory that can 
renormalize currents that are leading in v for a particular quantum number. To determine 
the NLL anomalous dimensions we compute the current correlators diagrams in Fig.[7|rather 
than two-loop vertex diagrams. The former method is more convenient since IR-divergences 
originating from the Coulomb phase cancel automatically and a distinction between IR- 
and UV-divergences becomes unnecessary |j| . In Fig. [7| four-quark interactions without 
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labels refer to the Coulomb potential and crossed circles to the currents (and their complex 
conjugate). Crossed squark lines refer to insertions of the NNLL order squark kinetic energy 
operator. Combinatorial factors and diagrams obtained by flipping the graphs left-to-right 
and up-to-down are understood. The results for the NLL order anomalous dimensions read 



d 



16vr 2 



+ ylP- + a\{mv) [3V«(i/) + 2Vg(i/)] . 



d 

-^Hc P {p)\ 



W/,a r^M 



48tt 2 



+ 



vi s) M 



+ a 2 (mu) 



Mi 



(31) 



The expression for the S-wave current agrees with the one for the S-wave current in the 
quark production case 0] (see also Refs. p, EH), up to the absence of the spin-dependent 

(T) 

terms and the fact that the coefficient V2 differs for the both cases. As the solution for 
Eqs. ()31|) we find 



In 



csQ) 



62 7ra s (m) (1 — z) + 63 ira s (m) \n(z) 



+ 64 7ra s (m) 



l-z 



l-13C A /(6/9o) 



+ 6 7ra s (m) 



z — 1 — w 1 ki(w) 



In 



cp\y) 

CP (1) 



where 



do 



dn 



di na s (m) (1 — z) + do na s {m) z—l — w 1 ln(w 



C F [C A C F {9C A - 100C F ) - P {26C 2 A + 19C F C A - 32C 2 ) 

gf [C A (9C A - 100C F ) - Q(3 (3C A - AC F ) 
2/? 2 (6/?o - 13CU) 

2AC 2 F (11C A - 3(3 )(5C A + 8C F ) 
13C A (6/3 - 13C A ) 2 

8C F (C A + C F )(C A + 2C F ) 
3/? 2 

C F (C A + 2C F ) 

m 

8C A C F (C A + AC F ) 
9(3 2 



(32) 



(33) 
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FIG. 8: NLL running of the normalized Wilson coefficient for S-wave and P-wave currents for 
m = 220 GeV (solid and dashed lines) and for m = 500 GeV (dot-dashed and dotted lines). 

In Fig. [Hj we have displayed the NLL running of the Wilson coefficients cs(v) and cp(v) 
normalized to their matching values. For the input parameters we have chosen two different 
values for the heavy scalar mass, m = 220 GeV and m = 500 GeV, and a s (mz) = 0.118, 
taking leading-logarithmic running for a s with rif — 5 active massless quark flavors and no 
active massless squarks (n s = 0). We find that the ^-variation for the S-wave coefficient 
is much stronger than for the P-wave coefficient. While the S-wave coefficient increases by 
about 6% for v ~ a s compared to the matching value at v — 1, the P-wave coefficient 
only increases by less than 3%. Furthermore, for S- and P-wave coefficients the maximum 
slightly decreases with the heavy scalar mass and also moves towards smaller values of v. 
The latter feature can be understood qualitatively from the fact that the average velocity 
(v) ~ a s (mQ s ) decreases with the heavy scalar mass m. 

V. CONCLUSION 

In this work we have presented the vNRQCD Lagrangian for a particle-antiparticle pair 
of colored scalars in the fundamental representation of QCD. The paper provides the needed 
ingredients for a NLL order description of the QCD effects of scalar-antiscalar production 
close to threshold in e + e~ and 77 collisions, and for scalar-antiscalar bound state quarks 
energies at NNLL order. 

The proposed Lagrangian is the spinless counterpart of the vNRQCD Lagrangian first 
formulated for quarks. To the order required in this work, the ultrasoft sector of the quark 
and scalar theories agrees, since to the order we are interested interactions with ultrasoft 
energy gluons are not sensitive to the spin. As a consequence the 1-loop ultrasoft renor- 
malization is identical in both theories. We also find that the soft coefficients in the scalar 
theory are in complete analogy with the corresponding ones in the quark theory. This is 
related to similarities in the structure of the HQET and HSET Lagrangians up to order 1/m 
due to reparametrization invariance. Differences arise in matching conditions that are not 
protected by symmetries and in the running induced by soft divergences since for scalars 
there are no mixing due to spin-dependent operators. 
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APPENDIX A: CONSTRUCTION OF THE TREE LEVEL HSET LAGRANGIAN 

The QCD Lagrangian for a colored scalar field with mass m reads 

£s Q CD = (D^YDy - m 2 <P*<p. (Al) 

To obtain the tree level heavy scalar effective Lagrangian from £qcd we need to take the 
limit m — > oo holding the heavy scalar's four- velocity t> M fixed (v 2 = 1). With this aim, the 
colored scalar field with momentum p = mv + k, where k is small, can be decomposed into 
two pieces, 

p —im v-x 

0(x) = — =^(^ + x ) , (A2) 
y 2m 

${x) = e imv - x P + <j>{x) , x(x) = e imv x P- <j>(x) , (A3) 



where 




(A4) 



project out the positive (particle) and negative (antiparticle) energy components of the field 
in the infinite mass limit. The explicit phase-factor introduced in Eq. flA3|) extracts the 
relativistic fluctuations of the particle field. This procedure is in complete analogy with the 
decomposition performed for a heavy quark field to write down the HQET Lagrangian, where 
the upper and lower components of the heavy quark field are projected out by P± = (l±^)/2. 
The \/2m factor has been introduced to achieve the non-relativistic normalization for the 
fields ip and x- 

The QCD Lagrangian with the decomposition above reads 



£sQCD = (if>* + X 



iv ■ D + 



{iv ■ Df D 2 



(V> + x) 



2m 2m 

= r(iv-D- ~X*(iv-D + 2m+ - X * - P ^X, (A5) 

where = — v^(v ■ D) is the component of the covariant derivative orthogonal to the 
velocity, i.e. v ■ D± = 0. In the second line of Eq. ()A5|) we have eliminated the term (iv ■ D) 2 
and the temporal derivatives in the crossed terms with the help of Eqs. (jA2HA4j) . (Doing the 
manipulations for ip and x separately also the equation of motion for the field in the full 
theory is required.) The field ip describes a massless degree of freedom, while x corresponds 
to fluctuations with twice the heavy scalar mass and can be integrated out of the theory at 
the energy scales we are interested in. At tree-level this can be done by solving the equation 
of motion for the field X'- 

2m + K ,.z> + |l) x = -|U, (A6) 
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which implies 



1 D\ , WD D\ \ D\ , 

2m + w-D + ^ 2m V 2m (2m) 2 J (2m) 2 



2m 

The expansion above is justified because derivatives acting on ip are of the order of the 
residual momentum of the heavy scalar, k = p — mv <C m. Inserting the relation (|A5|) in 
the QCD Lagrangian gives a local effective Lagrangian for heavy scalars as an expansion in 
1/m: 

* ( D\ D\ D\ D\ ( 1 \ \ 

HSET lp ^ 2 m g m 3 4 m 2 4 m 2 \ m 5 / J ^ ^ ^ 

The field redefinition 

(A9) 

can be used to eliminate the 0(l/m 4 ) time derivative term acting on ip and cast the tree 
level HSET Lagrangian in the form 



ip* ( iv 



2m 8m 3 I ^ ' ' 8m 2 ^ ' 4m 2 1 ( m 5 ) ) ^ 



•■HSET 

^* ■ £> - ^ + - g— ' y ±J +0[ A V>- (A10) 

1 2m 8m 3 y 32m 4 Vm 5 / 1 

In the heavy scalar rest frame, v = (1, 0, 0, 0), = (0, — D), the tree level HSET Lagrangian 
takes the form 



•■HSET 



„ D 2 D 4 {{D,E},D 2 } ^/1\\ /A x 



Note that at tree level the Darwin term [D ■ E]/(8m 2 ) is absent, and at order 1/m 3 only 
the kinetic D 4 /8m 3 operator exists. The full operator basis up to order 1/m 3 (relevant for 
HSET and HQET) can be found in Ref. 0. 
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